Self- Trapping of Polarons in the Rashba-Pekar Model 
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We performed quantum Monte Carlo study of the exciton-polaron model which features the self- 
trapping phenomenon when the coupling strength and/or particle momentum is varied. For the first 
time accurate data for energy, effective mass, the structure of the polaronic cloud, dispersion law, 
and spectral function are available throughout the crossover region. We observed that self-trapping 
can not be reduced to hybridization of two states with different lattice deformation, and that at 
least three states are involved in the crossover from light- to heavy-mass regimes. 
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PACS numbers: PACS numbers: 71.38.-k, 02.70.Ss, 71.38.Fp, 71.38.Ht 



Properties of particles strongly coupled to their envi- 
ronment are of importance in many fields of physics and 
are attracting constant attention given extreme diversity 
of what may be called a "particle," an "environment," 
and how they interact with each other. In most general 
terms, self-trapping (ST) means a dramatic transforma- 
tion of particle properties n when system parameters are 
slightly changed. LandaulJJ showed that the "trapped" 
(T) particle state with strong lattice deformation around 
it and the weakly perturbed "free" (F) particle state may 
have the same energy at some critical value of the cou- 
pling strength, a c . 

Of course, the resonance between F and T states is 
not infinitely sharp since the matrix element hybridiz- 
ing them is non-zero, i.e. ST phenomenon is a crossover, 
rather than a transition, and all polaron properties are 
analytic in a — see Ref. |] for an explicit proof. This the- 
orem makes the notion of ST rather vague since there is 
always some admixture of one state in another. More- 
over, it challenges the adopted opinion that only two, 
namely F and T, states are in competition. If there are 
more than two states within the energy scale of the hy- 
bridization matrix element then all of them are mixed 
and the F-T classification fails. In fact, the two-states 
assumption on which the current theory is based is not 
supported by experiments &»d rather complex spectra 
are usually observed insteadSo. 

According to the standard criterion™, ST takes place 
if there is a barrier Ub in the adiabatic potential be- 
tween the bare-particle and polaron states. It occurs, 
almost by definition, in the intermediate coupling regime 
where perturbation theory is not applicable. Hence, the 
existence of a barrier — if the very notion of the adiabatic 
potential is not ill defined — and the ST phenomenon can 
be addressed only by an exact method, because in the in- 
termediate coupling regime an analytic solution is hardly 
available, and even sophisticated variational treatments 
often give misleading resultaD. 

In this Letter, we consider a typical model in which 
particle couples to the environment of gapped disper- 
sionless optical phonons. For this model it is possible 
to define ST in a mathematically rigorous way and pro- 



ceed with its quantitative study. We show how various 
particle properties (energy, effective mass, dispersion law, 
and the structure of the polaronic cloud) change between 
weak- and strong-coupling limits, and provide detailed 
information about ST of polarons, which is not based on 
any approximations. Besides, we show that there are at 
least three states involved in mixing in the critical region 
and, thus commonly accepted concept of only F and T 
states mixing at a c appears to be oversimplified. In fact, 
we are not aware of any other numerical study testing 
how accurate are existing treatments of the ST problem. 

The Hamiltonian of the system consists of the free- 
particle term (we consider continuum three-dimensional 
case with dispersion relation e(k) = k 2 /2m ) 



h c = e ( k ) °k a k , 

k 

the Hamiltonian of the phonon bath 

q q 

and the standard density-displacement interactior 



ffc-ph = £ V(q) (&t - b _ 



(1) 



(2) 
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In Eqs. 



and b arc the particle and phonon anni- 



hilation operators in momentum space, correspondingly. 

Our study is based on the quantum Monte Carlo simu- 
lation of the polaron Green function in imaginary time at 
T = and.-a ubse quent analytic continuation to tr. 
frequencies0ll3llil. The method suggested in RefsBlMi 
is free from approximations and systematic errors. It is 
particularly suited for the study of ST problem where 
several (5-peaks are expected below the spectral contin- 
uum in the Lehman expansion 



S^(u;) = Y S{oj-E u (k)) \{v\a 



(4) 



Here {\v)} is a complete set of eigenstates of H in the 
momentum sector k, i.e. H |f (k)) = E v Qs) \v{k)). 
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Separating stable quasiparticle states (labeled by index 
i) from continuum, we rewrite Eq. (Q) as 

i J ^ 

(5) 

where Z^(Q) and are Z- factors and energies of 

stable states, and the continuum threshold is given by 
lu c = i?Q k_0 ^ + luq. Any state with E > lu c is unstable 
against single- (n — 1) or multi-phonon (n > 1) emission 
process nwn, where momentum p is selected 

only by the energy conservation law since phonons are 
dispersionless. 

Speaking rigorously, by self-trapping one understands 
the existence of such a region in the parameter space of H 
where more than one stable polaron states, differing by 
the degree of polarization of the lattice, coexist. This 
definition implies three critical points in the coupling 
constant (keeping other parameters fixed for simplicity), 
a c i(k) < a c (k) < a C 2(k). The ST "transition point" a c 
is understood as the point of avoided crossing between 
the two lowest polaron states. At this point the ground- 
state of the polaron is a hybrid of states with substan- 
tially different degrees of lattice polarization. The critical 
points a c i and a C 2 correspond to the appearance and dis- 
appearance of the extra stable state(s), respectively. [By 
definition, the energy difference AE^\a) between the 
ground and first stable excited state, which has its mini- 
mum at a c (k), ought to be less than to c — E^ .] Critical 
couplings introduced above are consistent with previous 
considerationsBfflj and have an advantage of being un- 
ambiguous, even when the minimal gap AE^ (a c (k)) is 
not smalM 

A typical system that is boliewed to feature ST is the 
so-called Rashba-Pekar modelld'E-a which describes Wan- 
nier exciton in the Is state interacting with optical vi- 
brations via electrostatic potentials 

V{q) = 7(q) { [1 + (La B q) 2 } 2 ~ [1 + (^a B q) 2 } 2 } ' 

(6) 

7 (q) = i (2V2«7r) 1/2 q' 1 . (7) 

Here a is the standard dimensionless coupling constant, 
cib is the Bohr radius, and £ e ,/i = TO e) ft/[2(m e + m^)] is 
given in terms of electron (m e ) and hole (m/,) masses, 
respectively. In this Letter, we focus on the [parameters 
corresponding to the curve number 2 of Ref.tJ, which is 
believed to describe ST in the strong-coupling regime. 
Setting the total bare mass of the exciton to = m e + mh, 
phonon frequency u>o, electric charge, and Plank con- 
stant to unity, one finds that m e = 0.065. The Bohr 
radius can be used to change the degree of adiabatic- 
ity in the model. Below we shall thoroughly consider an 
"almost adiabatic" case (i) with Ub /^>o — 2 (which is re- 
alised for as — 0.467) and outline some peculiar features 
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FIG. 1: The groundstate energy, average number of phonons, 
and effective mass as functions of a (points connected by solid 
lines). Relative statistic errors are less than 10~ 3 and 10~ 2 for 
the energy and (N) , respectively. The relative statistic errors 
for the mass are of order 10~ 2 for a < 18.5 and around 5 x 
10 -2 for larger coupling constants. Dashed lines show results 
of the perturbation theory while the dotted line corresponds 
to the strong-coupling limit. 

of the "nonadiabatic" situation (ii) with Ub/uq = 0.5 
(as = 0.934). The critical coupling constants, deter- 
mined within the approach of Ref.E^I, are then a^ d w 14.3 
and a^ d w 7.2, respectively. 

In Fig. 1 we show how the groundstate properties (k = 
0) depend on the coupling strength. The groundstate 
energy, the effective mass m* , and the average number 
of phonons in the polaronic cloud 

(7V) = (k = 0| £&t 6q | k = ), ( 8 ) 
q 

clearly indicate drastic changes around a c ~ 18.35. At 
this point the energy derivative changes very fast, and 
both (N) and m* ungergo step-wise increase visible even 
on the logarithmic plot for to*. In a narrow region be- 
tween a — 17.5 and a = 19 the effective mass increases 
by two orders of magnitude. The remarkable fact is that 
at a c the strong-coupling-, approach is still far from be- 
ing accurate (dotted lineE-3 in the upper panel of_Fig. [I]). 
Besides, the adiabatic critical constant o^ d w 21L3 differs 
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FIG. 2: Partial weights of n-phonon states in the polaron 
ground state (k = 0) at a — 18 (circles), a = 18.35 (squares), 
and a = 19 (dimonds). Statistic errorbars of order 3 x 10~ 3 
are less than the symbol size. 



significantly from our value a c ~ 18.35. For the "nonadi- 
batic" case (ii) the behavior of the groundstate proper- 
ties is qualitatively the same, but quantitative deviations 
from the strong-coupling limit are larger. 

Next, we study how the phonon cloud evolves through- 
out the ST critical region. Partial n-phonon contribu- 
tions to the polaron ground state Z^(n) are the prob- 
abilities of finding exactly n phonons in the cloud, and 
the average number of phonons introduced earlier, is just 
(N) = Xi n nZ(n). Figure ^] shows Z^ =0 (n) distributions 
at a = 18 (below the crossover region), a = a c — 18.35, 
and a = 19 (trapped state). We see that the distribu- 
tion at a c has two peaks and is half-way between the 
two limiting cases. However, in the "nonadiabtic" case 
(ii) the structure with two maxima in Z{n) is missing. 
Therefore, the peculiar behavior presented in Fig. ^ is a 
general feature of the ST phenomenon whereas the two- 
peak structure of the phonon distribution is specific for 
the adiabatic limit. 

The spectral function around the critical pointQ re- 
veals up to three stable excited states below the contin- 
uum threshold (see examples of the Lehman function in 
Fig. ||) . We observe in Fig. ^ that three polaronic states 
(in the energy range comparable with the hybridization 
strength) participate in the ST crossover. We underline 
that all three states have large Z-factors (> 0.1). There- 
fore, for the given set of parameters more than two states 
are mixed at the crossover point and the standard pic- 
ture of F-T hybridization at the tip of the ST crossover 
fails. One can speculate that extra stable states in the 
gap, which standard theory puts into the spectral con- 
tinuum, are due to excited levels of highly nonlinear ST 
potential in the resonating region. However, this inter- 
pretation is essentially qualitative since the concept of 
adiabatic potential breaks down in the crossover region. 
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FIG. 3: The Lehman spectral function 5'' k=0 - ) at coupling 
constants a = 18.35 (upper panel) and a — 18.75 (lower 
panel) . 
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FIG. 4: Energies of the ground (circles) and excited sta- 
ble states (squares, diamonds, and triangles) vs interaction 
constant. The dashed line is the threshold of incoherent con- 
tinuum. Typical errorbars for the first, second, and third 
excited states are 10 -2 , 3 x 10~ 2 , and 4 x 10~ 2 , respectively. 



So far we have considered ST crossover at zero momen- 
tum, i.e., for the ground state. However, same considera- 
tions apply to finite momentum states, as long as k < k c 
where k c is defined as the point whcxoJi , Q k ' = uo c and the 
polaron spectrum has an end pointBEj. For a < a c the 
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FIG. 5: The wave- vector dependence of energy and aver- 
age number of phonons for a — 17.75 (circles connected by 
solid lines). The errorbars are 3 x 10" 3 and 10~ 2 for energy 
and (N), respectively. The dashed curve is the effective mass 
approximation _B (k) = E + k 2 /2m* with E (a = 17.75) = 
-3.7946 and m*(a = 17.75) = 2.258 obtained from direct 
Monte Carlo estimators. The dotted curve is the parabolic 
dispersion law fitted to the last four points in the energy plot 
with parameters E\(a = 17.75) = —3.5273 and m\ = 195. 
The open rectangle is the energy obtained from spectral anal- 
ysis for the first excited state. 



ground state is characterized by small lattice distortion 
and light effective mass and thus the dispersion law asso- 
ciated with this state rises steep with k. On another 
hand, the dispersion curve of the excited heavy-mass 
state is nearly flat which means that light- and heavy- 
branches have to intersect at some wave vector. In Fig. ^ 
we plot the groundstate dispersion law and the average 
number of phonons in the ground state for a = 17.75 
which we interpret as the level crossing picture (see the 
figure caption). Notice the agreement between the spec- 
tral analysis at k = and the parabolic fit of the heavy 
mass branch. 

By all accounts, the existence of more than one stable 
polaron state is a highly nontrivial qualitative property of 
the model whether the hybridization gap between these 
states is small, or not. Moreover, the ST crossover is not 
necessarily limited to hybridization of only two states. 
Recent studies of the Holstein polaron in ID strongly sup- 
port the universality of the latter statement since more 
than two stable states were found there as wellliJ. 

We found that the dependence of Eq and (N) (but not 
the structure of the cloud) on coupling can be used as an 
indirect indication of the self trapping phenomenon with 
a well-defined "transition point" (the point of minimal 
gap between the ground and first stable excited state). 
It was considered previously as a, .theorem, that ST may 
occur only in dimensions d > 2du; however, the defini- 
tion of what has to be counted as a ST transition was not 
given in quantitative terms. We believe that it is more 
appropriate to use less "radical" , but unambiguous, defi- 
nition adopted in our paper. Recently, the second stable 
polaron state was found to exist for the Holstein polaron 
in a one-dimensional latticaia and in infinite dimension 
appr oximat ionlI3 . 
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